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ABSTRACT 

Rapidly accreting massive protostars undergo a phase of deuterium shell burning during pre-main 
sequence evolution that causes them to swell to tenths of an AU in radius. During this phase, those 
with close binary companions will overflow their Roche lobes and begin transferring mass. Since 
massive stars frequently have companions at distances well under 1 AU, this process may affect the 
early evolution of a substantial fraction of massive stars. We use a simple protostellar evolution model 
to determine the range in accretion rates, mass ratios, and orbital separations for which mass transfer 
will occur, and we compute approximately the stability and flnal outcome of the transfer process. We 
discuss how mass transfer affects the demographics of massive binaries, and show that it provides a 
natural explanation for the heretofore unexplained population of massive "twins" , high mass binaries 
with mass ratios very close to unity. 

Subject headings: accretion, accretion disks — binaries: close — binaries: spectroscopic — stars: 
evolution — stars: formation — stars: pre-main sequence 



1. INTRODUCTION 

Massive stars form in regions of high pressure and den- 
sity that produce high accretion rates. Simple order- 
of-magnitude estimates for the properties of observed 
massive protostellar cores suggest that accretion rates 
of ~ 10-"* - 10-2 Mq yr-i are typical (Krumholz 2006'), 
and both detailed analytic m odels ([M cKcc & Tan 2003) 
and numerical simula tions (jBanerjee fc Pudritz. ,2007i : 
iKrumholz et al.l |2007[ ) produce accretion rates in this 
range. Observations of protostellar outflows with mass 
fluxes up t o ~ 10~2 MfD yr~^ f rom luminous embedded 
protostars (jHenning et al.ll2000D also suggest high accre- 
tion rates. 

The combination of rapid accretion and high mass pro- 
duces protostars with very large radii, for two reasons. 
First, Stabler (1988) shows that prior to the onset of 
deuterium burning, the radius of an accreting protostar 
is determined primarily by the specific entropy of the 
gas after it passes through the accretion shock, radiates 
and settles onto the surface, and is eventually buried in 
enough optical depth to prevent it from radiating fur- 
ther. A high accretion rate reduces the amount of time 
an accreted gas element has to radiate before it is buried, 
and thus prod uces higher specific entropy an d larger ra- 
dius. Second, iPalla fc Stahle^ (|1991l . Il992t l show that 
after deuterium burning begins, massive protostars pass 
through a period of shell burning. This occurs because, 
as contraction raises the core temperature, the opacity 
decreases to the point where the core becomes convec- 
tively stable. Accreting deuterium cannot penetrate the 
radiative core and, as a result, deuterium burning occurs 
in a shell around it. In a process analogous to that in a 
red giant, this produces rapid expansion of the envelope 
above the burning layer. Together, these two effects can 
produce radii of several tenths of an AU during pre-main 
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sequence evolution of rapidly accreting stars. 

Large radii create the possibility for mass trans- 
fer in close binaries. Most massive stars appear 
to be in m ultiple s ysterns (e.g. [Preibisch et al.l 120011 : 
IShatskv fc Tokovini^[200l |Lada| 120061) . While the semi- 
major axis distribution of massive stars is not well- 
determined, due to low statistics and complex selec- 
tion biases, there appears to be a significant popu- 
lation of very close binaries. The WR20a system, 
the most mass ive binary known, has a separa.tion o f 
only 0.25 AU (jBonanos et al.l 120041: iRauw et al.ll2005l ). 
The semi-major axis of the massive detached eclips- 
ing binary D33J013346.2 -H3 04439.9 in M33 is 0.22 
AU (iBonanos et all l2006l). iHarries et all (|2003l) and 
iHilditch et al.l ( 20051 ) report a sample of 50 OB star 
eclipsing binaries with periods of 5 days or less taken 
from the Optical Gravitational Lens Experiment sur- 
vey (OGLE; ildalski, et al. 1998) of the SMC. This pe- 
riod limit corresponds to a semi-major axis of 0.2QM\I^q 
AU, where Mioo is the total mass in units of 100 Mq. 
Roughly half the systems are detached binaries, indicat- 
ing that they have not yet undergone any post-main se- 
quence mass transfer. It is therefore likely that these 
systems formed at close to their current separations, or 
possibly even closer, since conservation of energy dictates 
that mass loss from winds during main sequence evolu- 
tion widens the orbits of massive binaries. This means 
that systems such as WR20a, DD33J013346.2-H304439.9, 
and the SMC binares almost certainly experienced a 
phase of mass transfer during their pre-main sequence 
evolution. 

To study how this process will affect massive binaries, 
we proceed as follows. We assume that a "seed" close 
protobinary system forms while the system is still em- 
bedded and accreting. Such a binary may form via one 
of two possible mechanisms that have been proposed: di- 
rect fragmentation of a massive molecular cloud core and 
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around a massive protostar and subs equent migration of 
a ste l lar-mass fragment inw ard (e.g. iKratter fc Matzneil 
120061: iKrumholz et al.ll2007[ ). Although these two mecha- 
nisms arise in the context of specific star formation mod- 
els, we emphasize that seed binaries of this sort are ex- 
pected to form early in massive protostars' lives in the 
context of any model where massive protostars form by 
accretion (e.g. Mc Kee fc TmI 120031: iKeto fc WoodI 120061 : 
Krumholz ct al. 20Q3)- Moreover, recent observations by 
Apai et al.l 12003) show that massive stars have a high 



intrinsic binary fraction even while they are still embed- 
ded, providing observational support for this hypothesis. 
It is worth noting, however, that stable binaries may not 
form until much later in massive pr otostellar evolution if 
massive stars form by collisions (e.g. Bon nell et al.|[l998l ). 
since binary orbits will suffer continual disruption as long 
as close encounters leading to mergers are occuring. If 
massive stars form primarily by collisions and not by ac- 
cretion, the calculations we present in this paper will not 
apply. 

Given a seed binary, in this paper we calculate how pre- 
main sequence stellar evolution and mass transfer will 
modify its properties. In § [21 we determine the range 
of parameters for which mass transfer can occur, and in 
§ [3] we discuss the outcome of mass transfer when it does 
occur. We discuss the implications of this process for the 
massive binary population, and related points, in §21 and 
we summarize our conclusions in § [5l 

2. THE PROPERTIES OF MASS TRANSFER BINARIES 

The problem of binary pre-stellar binary evolution de- 
pends on both the accretion histories of the stars and on 
the properties of their orbit. Even if we limit ourselves to 
protostars accreting at constant rates and moving in cir- 
cular orbits of constant size, there are four parameters: 
the accretion rates onto each of the two stars, the offset 
in time between formation of the two stars, and the semi- 
major axis of the orbit. Since fully exploring this param- 
eter space would be time consuming and not particularly 
informative, we further simplify it by assuming that the 
two binary components form coevally or nearly so, re- 
ducing the space to three dimensions. Our goal is simply 
to map out for what ranges of accretion rate and orbital 
separation mass transfer is a possibility. Our assump- 
tion that the two components of the binary are roughly 
coeval and that the larger component has a higher time- 
averaged accretion rate (necessarily the case for coeval, 
constant accretion rate systems) is sup ported by recent 
simul ations of massive star formation (jKrumholz et all 
[2001 . 

The first step in exploring this parameter space is to 
determine the protostellar mass-radius relation as a func- 
tion of accretion rate, which we do in § 12.11 Using these 
tracks, in § 12.21 we determine the minimum semi-major 
axis required for Roche lobe overflow (RLOF) to occur 
in a binary system consisting of two stars of mass Mi 
and M2, Ml > M2, with secondary-to-primary mass ra- 
tio q = M2/M1 < 1, orbiting with semi-major axis a. 
The total system mass is M, and the binary accretes at 
a rate Mace- Under our assumption of coeval formation, 
the ratio of accretion rates onto the two components is 
-^•f2,acc/Afi,acc = M2/M1 — Q, whcrc the subscript "acc" 
indicates the change in mass due to accretion into the 
system, rather than transfer between the two stars. 



2.1. The Protostellar Mass-Radius Relation 

For a given accretion rate onto a star we construct a 
track of radius versus mass using a simple protostellar 
model. To allow quick exploration of a lar ge range of pa- 
ramet ers, we use the one-zone model of iMcKee fc TanI 
(|200l to generate our tracks. This model has been 
calibrated again s t the detailed numerical calculations of 
[Palla & Stahler (1992), and agrees to ~ 10%. We re- 
fer readers to McKee fc TanI for a detailed description 
of the model, but here we summarize its most impor- 
tant features. In this model a protostar is assumed to 
be a polytropic sphere with a specified accretion rate as 
a function of time. At every time step of model evo- 
lution one uses conservation of energy to determine the 
new radius, including the effects of energy lost in disso- 
ciating and ionizing incoming gas, energy radiated away, 
and energy gained by Kelvin-Helmoltz contraction, deu- 
terium burning, and the gravitational potential energy 
of accreting gas. 

An accreting massive protostar passes through four dis- 
tinct phases of evolution before reaching the zero-age 
main sequence (ZAMS). When the star first forms, it 
evolves passively without any nuclear burning. When its 
core becomes hot enough (~ 10^ K), the second phase 
begins: deuterium ignites, the star becomes fully con- 
vective and begins burning its accumulated deuterium 
reserve, and contraction of the core temporarily halts. 
Once the star exhausts its deuterium reserve, it enters 
a third phase in which it continues burning deuterium 
at the rate it is brought in by new accretion. In this 
phase the core resumes quasi-static contraction. As the 
core temperature continues to rise, its opacity drops, and 
eventually a layer forms that is stable against convection, 
initiating the fourth phase. Newly accreted deuterium 
cannot pass through the radiative layer to reach the core, 
and without an influx of new deuterium the core rapidly 
exhausts its supply and becomes convectively stable as 
well. Deuterium begins burning in a shell above the ra- 
diative zone, driving a rapid expansion of the star. As the 
core continues contracting, the radiative zone increases 
in size and the star resumes contraction. Eventually the 
core becomes hot enough to ignite hydrogen, and at that 
point the star joins the ZAMS. 

Figure [1] shows some sample tracks of radius versus 
mass computed using this model. The sharp increase in 
radius that each model shows is the result of the start 
of deuterium shell burning. As the plots show, for the 
high accretion rates expected in massive star-forming re- 
gions, the radius can reach several tenths of an AU. We 
compute protostellar models for 600 values of Mace, uni- 
formly spaced in logarithm in the range 10~^ — 10~^ M© 
yr^^, accreting up to a maximum mass of 100 Mq. By 
this mass, stars have joined the ZAMS regardless of their 
accretion rate. 

2.2. Roche Lobe Overflow 

Having determined the protostellar mass-radius rela- 
tion, we can now determine for what semi-major axes 
a a binary of mass ratio q (secondary-to-primary) will 
experienc e RLOF. The r adius of the Roche lobe around 
star 1 is ([EggletonlllQSl) 

_ 0.49 

"^'^^ 0.6 + 52/3 iog(l + g-i/3)' 
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Fig. 
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— Radius versus mass for constant accretion rates of 
2 Mp 



Ipi yr ^ (solid lines), computed using the simple one- 
zone IMcKee fc Taii ^003) model. For comparis on we also show 
the results of the detailed numerical calculation of JPalla fc Stahleil 
|[1991 . 1992) for accretion rates of 10"^ Mq yr"! and 10"" Mq 
yr~^ {dot-dashed line). The sharp rise in _R that all models show 
results from the start of shell deuterium burning. 

the Roche lobe radius around star 2 is given by an anal- 
ogous expression with q replaced by q^^. For a given 
value of q and Mi^acc, the accretion rate onto the more 
massive star, it is straightforward to use the mass-radius 
tracks to compute whether RLOF ever occurs for a given 
value of a, and, if so, to determine the mass and state of 
each star when it does. 

Figure [5] summarizes the results for 3 different values 
of q. We find that in rapidly accreting systems, RLOF 
will occur when the primary star is Mi ^ 10 Mq and the 
semi-major axis a is a few tenths of an AU or less; a 
rough analytic estimate for the critical value of a below 
which RLOF occurs is 



0.2 My/(1 



AU, 



(2) 



where M_4 = Mi_acc/10~** Mq yr~^ Thus, in the con- 
text of the model proposed here, observed systems such 
as WR20a, D33J013346. 2+304439.9, and the eclipsing 
binaries in the SMC, with masses > 20 Mq and separa- 
tions < 0.25 AU, are well within the expected overflow 
range for reasonable accretion rates and initial mass ra- 
tios. We further find that it is always the more massive 
star that overflows its Roche lobe, rather than the less 
massive one. This is because, despite the larger Roche 
radius around the more massive star, its higher accretion 
rate leads this star to have a larger radius before the on- 
set of shell burning, and to start deuterium shell burning 
sooner. Finally, note that much of the parameter space 
allowed for transfer, the region between the dashed and 
solid lines in Figure O is due to the shell burning phase, 
which roughly doubles the radius of the star and there- 
fore doubles the minimum value of a for which transfer 
can occur. 

If we were to relax our assumption that the stars 
are coeval, our conclusion would be strengthened, be- 
cause simulations indicate that in massive binary sys- 
tems the more massive companion generally forms earlier 
(jKrumholz et al.l 120071 ). Since the mass at which deu- 
terium shell burning starts increases with accretion rate 
(see Fig. [Ij , and the protostellar radius at any time in- 
creases with both mass and accretion rate, this general 
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Fig. 2. — Minimum a for RLOF {solid line) as a function of 
accretion rate onto the more massive star Mi ^cc. The value of q 
is indicated in each panel. Below the RLOF line, in the region of 
parameter space where overflow occurs, we show the mass Afi (in 
solar units) of the donor star at overflow {shaded regions). Binaries 
with a below the dashed line experience overflow before the primary 
starts deuterium shell burning, while those with larger values of a 
experience overflow after the onset of shell burning. 

result should be robust. 

3. THE OUTCOME OF ROCHE LOBE OVERFLOW 

3.1. Time Scales 

Before we attempt to calculate the outcome of RLOF, 
it is helpful to review some relevant timescales for the 
problem. The shortest is the orbital period, P = 

1^/2 1/2 

3.7aQ Mj^Q days, where ao.i is the semi-major axis 
in units of 0.1 AU and Miq is the total system mass in 
units of 10 Mq. This is the time scale on which mass 
lost from the primary star will reach the vicinity of the 
secondary; actual accretion may take longer if the gas 
has to be processed through a disk. Next is the dynami- 
cal, sound-crossing time of the donor star. The sound 
speed varies greatly from the stellar core to the sur- 
face, so to be conservative and obtain the longest pos- 
sible time scale we evaluate this using the surface sound 

speed, which gives Tdyn = -R*i/cs,surf ^ 31i?5oT'4 
days, where i?5o is the stellar radius in units of 50 Rq 
and T4 is the stellar surface temperature in units of 10* 
K. The timescale Tdyn describes the time required for 
the star to adjust mechanically as it loses mass. Finally, 
there is the Kelvin-Helmholtz time of the donor star, 
TKH GMl/{R^iLi) - 62007\ffoi?5o^.L4 1 yr, where L4 
is the star's luminosity in units of 10'* Lq. This describes 
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the time required for the star to adjust thermaUy to mass 
foss. 

The most important point to take from this calculation 
is that tkh is by far the largest timescale in the problem, 
so that stars will be unable to adjust their thermal (as 
opposed to mechanical) structure to mass transfer that 
occurs on orbital or dynamical timescales. Thus, we may 
approximate stars as behaving adiabatically during mass 
transfer. A corollary of this is that, as gas in the donor 
star expands adiabatically in response to mass loss, deu- 
terium burning will slow dramatically, since the burning 
rate is extremely temperature sensitive. Thus, once the 
star overflows its Roche lobe, its envelope will not have 
its specific entropy altered by further nuclear burning on 
a dynamical time scale. 

3.2. Stability of Mass Transfer 

We first address the question of whether the RLOF 
leads to stable or unstable mass transfer. Transfer is 
stable if adiabatic mass loss shrinks the star at a rate 
faster than the Roche lobe around it shrinks, and is un- 
stable otherwise. If the transfer is unstable, the donor 
star can change its mass by order unity on a time scale 
Tdyn, and transfer will stop only when a new mechanical 
equilibrium is established. 

To check stability, we must evaluate how mass trans- 
fer changes both the Roche radius and the stellar radius. 
The former is easy to compute. If no mass is lost from 
the system, then conservation of angular momentum de- 
mands that the radius of the system shrink at a rate 



(3) 



a Ml ( 1 
- = -2— 1 

a Ml \ go, 

where qo is the secondary-to-primary mass ratio at the 
time when mass transfer begins and Mi is the rate of 
change of the primary's mass due to mass transfer (as 
opposed to accretion onto the system). Since Afi < 
and qo < 1, this means that d < 0; the semi-major axis 
of the orbit shrinks in response to mass transfer. From 
equation dT]), the Roche radius around the primary star 
varies approximately as 
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The expression in square brackets is always positive for 
qo < 1, so the Roche radius around the donor star shrinks 
as well. 

For the range of parameter space we explore, a star un- 
dergoing RLOF does so before reaching the ZAMS, and 
at a mass below 30 Mq. Thus, the star is gas pressure- 
dominated, at least in its outer layers. If RLOF occurs 
before the onset of deuterium shell burning, the star is 
convective throughout, and thus is well-described as an 
isentropic n = 3/2 polytrope. Since a star with this 
structure expands in response to mass loss at constant 
entropy (see below), we learn immediately that RLOF 
will be unstable in this case, because the star expands as 
the Roche radius contracts. Thus, systems whose param- 
eters fall below the dashed lines in Figure [5] are always 
unstable. 



If RLOF occurs after deuterium shell burning starts, 
the star consists of a lower specific entropy inner radiative 
zone and a significantly larger, high specific entropy en- 
velope. Shell burning eventually raises the temperature 
in the envelop e to the point w here it becomes radiative 
as well (Palla fc Stahleilll991f l. but at the onset of shell 
burning it is still convective. RLOF occurs at the start of 
shell burning, since this is the point of maximum radius, 
so the envelope will also be an isentropic n — 3/2 poly- 
trope, with a different specific entropy than the interior 
radiative zone. 

Foll owing the analysis of iPaczvhski fc Sienkiewicd 
()1972f ) for a red giant, we can approximate this configu- 
ration in terms of a "centrally condensed" polytrope, a 
polytropic sphere with a point mass in its center. This 
is a less accurate approximation than it is for a red gi- 
ant, since deuterium shell burning in a protostar only 
expands the envelope by a factor of slightly more than 
2, as opposed to more than an order of magnitude for a 
red giant. However, we can regard the pure polytropic 
model of a star before shell burning and the centrally 
condensed polytropic model for a star during shell burn- 
ing as limiting cases. 

Centrally con densed n = 3/2 polytropes obe y a mass- 
radius relation (jPaczvnski fc Sienkiewicj|1972r ) 



(5) 



where i^T is a constant that depends on the specific en- 
tropy in the convective envelope and _E is a dimensionless 
number that depends only on ^, the ratio of the core mass 
to the total mass. Since the star behaves adiabatically as 
mass is removed, K remains constant. If the mass of the 
radiative core is also unchanged, a reasonable approxi- 
mation when transfer begins, then the radius of the star 
varies as 



R*i 



Ml 



1 



Co 



mo) 



dE 



(6) 



where ^o is the core mass ratio at the start of mass trans- 
fer. The star expands for < 0.22, and shrinks for 
Co > 0.22. 

Since Rri = R*i at the start of mass loss, mass transfer 
is unstable unless < -R^i- This condition is met only 
if 



(7) 




(1 - go) 



B y evaluating E and dE/d^ nu merically (e.g. table 1 
of iPaczvnski fc Sienkiewic3ll972f ). this condition lets us 
identify for each ^o what values of the initial mass ratio 
qo produce stability. Figure [3] shows that for ^o < 0.22 
no stability is possible because the star expands rather 
than contracting in response to mass loss. 

Stars without radiative cores have ^ = 0, so if RLOF 
occurs before the onset of shell burning then ^o = and 
the transfer is unstable. Shell burning stars will have 
Co > 0, and they may be stable if Co and go are large 
enough. The value of Co is determined by the size of the 
radiative core when it first appears, which is in turn set 
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Fig. 3. — Minimum initial mass ratio go required to make mass 
transfer stable for a given initial core mass fraction ■ The dotted 
line at go = 0.22 indicates the minimum go below which mass 
transfer is always unstable regardless of the value of qo . 



by the radius of the sheU where the maximum luminosity 
that can be transported by radiation, Li-i,d{m), first rises 
to the point where it is equal to the luminosity passing 
through that shell, L{m). Here, m is the mass enclosed 
within a given shell, and varies from m = at the center 
of the star to m = Mi at its surf ace. 

From their numerical models, iPalla fc Stahled (|1992l ) 
find that Lra,d{m) = Lira) is first satisfied at m/Mi = 
Co = 0.20, 0.2f, and 0.41 for accretion rates of 10-^ 
3 X 10~^, and lO""* Mq yr"^ We would hke to extrap- 
olate to higher accretion rates, but that is quite difficult 
because varies so non-linearly with Afi^acc- Instead, 
we make some general observations to help understand 
the likely value of ^q- The non-linearity in the varia- 
tion of Co with Ml acc result s from the complicated shape 
of Lrad(m) (see figure 5 of IPalla fc Stahled UMI) . This 
shape admits two possible solutions to irad(w) = i(m), 
i.e. two mass shells where a convectively stable zone could 
appear, one at smaller m and one at larger m. The rea- 
son Co is nearly unchanged between Mi^acc — 10^^ and 
3 X 10~^ Mq yr~^ is that L^s.d{m) = L{m) is first sat- 
isfied at the same solution point in the two models; Co 
jumps between Mi, acc = 3 x 10~^ and 10""* Mq yr~^ 
because the point where Lrad(w) and L{m) first become 
equal moves from the inner to the outer possible solu- 
tion. Since the shape of the irad('Ti) curve does not sug- 
gest that there are other, larger mass solutions where 
Liadi'm) = L{m) might first occur, it seems unlikely that 
Co will be much larger than 0.41. However, we caution 
that this is a tentative conclusion, and that one cannot 
confidently estimate Co without detailed numerical mod- 
eling. In particular, t he va lue of Co hi the models of 
IPalla fc Stahleil (|1991l . I1992D appears to be sensitive to 
their choice of boundary condition. 

If our tentative conclusion holds and Co ^ 0.5 regardless 
of the accretion rate, then transfer will be unstable unless 

go ^ 0.8. 

3.3. The Fate of Stable Mass Transfer Binaries 

Massive binary systems that initiate stable mass trans- 
fer lie in the region of parameter space above the line in 
Figure [21 Stable transfer is possible only if it starts be- 
cause the more massive star has initiated shell burning, 
and possesses a radiative core. While transfer reduces the 
temperature in the envelope via adiabatic expansion and 
thus halts shell burning immediately, this seems likely to 



produce only a temporary reprieve. As the primary con- 
tinues to contract toward the main sequence on a Kelvin- 
Helmoltz timescale, it cannot avoid re-heating and again 
igniting deuterium in its envelope, producing further ex- 
pansion. Thus, there may be multiple episodes of shell 
burning, expansion, and stable mass transfer separated 
by time intervals of order tkh , each leading to the trans- 
fer of a relatively small amount of mass but over time 
accumulating to push the system considerably toward 
q = 1. Determining the detailed evolution of a system 
undergoing this process requires the ability to follow the 
internal structure of the donor star as it repeatedly loses 
small amounts of gas from its envelope, cools adiabati- 
cally, and then heats up again. This is beyond the ca- 
pabilities of our simple protostellar evolution model, and 
since our best estimate for Co indicates the unstable mass 
transfer is likely to be a more common occurrence than 
stable transfer, we do not pursue the problem further in 
this paper. 

3.4. Termination of Unstable Mass Transfer 

We have found that unless the mass ratio is already 
near unity and the donor star is shell burning, mass 
transfer is unstable and proceeds on a dynamical time 
scale. This makes it difficult to determine the precise 
outcome without a protostellar model that, unlike ours, 
is capable of following the accretion of very high en- 
tropy gas on this time scale. Numerical simulations of 
the transfer itself would also be helpful, since in the case 
of runaway transfer it is possible that mass lost from 
the primary may go into a circumbinary disk or a wind 
rather than accreting onto the secondary. Nonetheless, 
we can make some simple calculations to estimate how 
and whether mass transfer will cease. 

Mass transfer will continue until either the primary 
shrinks within its Roche lobe, or until the secondary over- 
flows its Roche lobe too. We discuss the latter outcome 
in § 13.51 We can check whether the former is a possibility 
by asking whether there exists a new mass ratio g, larger 
than the original mass ratio goj such that the primary 
and secondary are both within their Roche lobes. Since 
such an equilibrium must be established on the dynam- 
ical time scale, the donor star must be able to reach it 
adiabatically. To determine whether such an equilibrium 
exists, we first examine how the Roche radius and stellar 
radius for each star change with q in § 13.4. 1[ and then 
we search for equilibria in § 13.4.21 

3.4.1. Roche Radius and Stellar Radius 

If the mass transfer conserves both total system mass 
and angular momentum, then the ratio of the semi-major 
axis a to its value oo before transfer starts is 



Go 



1 



go 



(8) 



when the mass ratio reaches g. From equation ([T]), the 
new Roche radius of the primary is related to the original 
one Rri^o by 



Rrl 
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t1,0 



1 + go 
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2/3, /, , -l/3x 

go log(l + go ) 
0.6-Hg2/3 1og(l-Fg-i/3) 
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The ratio of current to initial Roche radius for the second 
star, Rr2/ Rr2,o, is given by an equivalent expression with 
q and go replaced with and q^^^, respectively. Note 
that once mass transfer on a dynamical timescale begins, 
the two stars cannot remain exactly synchronous, nor 
can they be treated as point masses, and thus equations 
([8]) and ([9]) are approximate (J. Goodman, 2007, private 
communication) . 

Similarly, from the mass-radius relation ([5]) , the initial 
radius of star 1 and its radius once the mass ratio reaches 
q are related by 



R 



1 



1 



90 



1/3 



E 



1+90 



^0 



2/3 



(10) 



This relation will break down if ^ = ^o(l + 9)/(l + '?o) be- 
comes too large, because formally a centrally condensed 
polytrope goes to zero radius as ^ — > 1. In reality, the 
ra diative core has a finit e radiu s, which the simulations 
of lPalla fc Stahleil (Il99ll . Hool) show is ~ 1/2 of a dex 
smaller than the total stellar radius at the onset of shell 
burning. However, we find below that the radius of the 
primary star can never become this small before the sec- 
ondary overflows its Roche lobe, so equation (jTU]) is an 
adequate approximation for our purposes. 

Now consider the secondary star. Mass transfered from 
the primary to the secondary will not have time to cool 
and will therefore retain its specific entropy; it may even 
gain entropy due to an accretion shock as it falls onto 
the secondary. Since this high entropy material is being 
transferred to the smaller potential well of the second 
star, it will rapidly expand the secondary's radius, and 
eventually cause it to overflow its Roche lobe as well. 

We can obtain a rough lower bound on the secondary's 
radius by treating it as a centrally condensed polytrope 
as well. In this approximation, we take the secondary's 
radius at the onset of transfer to be negligible in com- 
parison to its radius after mass transfer starts, and thus 
treat the secondary at the start of transfer as a point 
mass. After an amount of mass AM has been trans- 
fered, the secondary has mass M2 = M2.0 + AM and 
envelope-to-core mass ratio ^2 = AM/M2,o- Since the 
specific entropy of the gas that is transferred is the same 
as or larger than it was before the transfer, the constant 
K for the secondary is at a minimum the same as for the 
primary at the onset of transfer, 

K = i?,i,oA/i/o'i^(eo)"'/' - Rri^oM^i^Ei^o 



-2/3 



(11) 



The radius of the secondary is therefore bounded below 

by 



i?*2 > R*1.0 



1 



1 1/3 I E 



(i + go)g 



90 ( J±2.^ 

q V i+go J 



mo) 



(12) 



at the point when the mass ratio is q. 

3.4.2. Existence of Mechanical Equilibria 

Having determined how the Roche radius and stellar 
radius for each star vary as the stars exchange mass and q 
changes, we are now in a position to determine whether 
there is a state the system can reach adiabatically in 
which neither star is overflowing its Roche lobe. Such an 
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Fig. 4. — Ratio of Roche radius Rr to stellar radius Rt for the 
primary star (solid lines) and the secondary star (dashed lines) 
for a system with qo = 0.5 initially. We show curves for = 
0, 0.25, 0.5, and 0.75, from lowest to highest. A system stabilizes 
if there is a value of q for which Rr/R* > 1 for both stars. The 
asterisks mark an example of such a stable point: a system with 
go = 0.5, Jo = 0.75 stabilizes at g = 0.54. 




Fig. 5. — Value of q for which mass transfer terminates, as a 
function of Jo ^'nd 90- Outside the shaded region, mass transfer 
is either always stable [top right region), or there is no value of q 
for which it will terminate (bottom left region). Inside the shaded 
region, the grayscale shows the value of q for which mass transfer 
terminates. The lightest region corresponds to g > 0.8, the next 
darkest region to 0.6 < g < 0.8, the next to 0.4 < g < 0.6, and so 
forth. Althou gh Jo is uncertain, in i|3^2l we argue on the basis of 
the models of lPalla fc Stahleri 1 119911 . Il992') that it is likely to be 
< 0.5. 

equilibrium exists if and only if there is a value q > q^ 
such that R.^,l < Rj-i and R^,2 < Rr2- We also require 
q < 1. Since the envelope of the primary has higher 
entropy than that of the secondary, the secondary will 
certainly overflow its Roche lobe by the time g = 1 is 
reached, even though our lower limit in equation (|12p 
does not reflect this. We illustrate an example of the 
problem graphically in Figure [4l which shows the ratio of 
Roche radius to stellar radius for both stars as a function 
of q for go — 0.5 and various values of ^o- 

As the figure shows, for a given value of qo and 1^0 there 
may or may not be a value of q satisfying the condition 
that both stars are within their Roche radii. If there 
is such a q, then mass transfer will halt once the system 
reaches that value. We summarize these results in Figure 
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[5l which shows the value of q for which a system stabiHzes 
as a function of qo and ^o- In the upper right white region 
of the pfot, systems are initially stable and mass transfer 
does not occur on a dynamical time scale (see ij 13. 3p . In 
the lower left region, systems are always unstable, and 
there is no value of q for which the the donor star can 
shrink within its Roche lobe. Systems whose parame- 
ters fall within the shaded region in between are initially 
unstable, but can stabilize and stop transferring mass 
for some value of q. Generically, we find that systems 
are unable to reach equilibrium unless they fall within 
a relatively small region of parameter space. Note that 
even this certainly overestimates the size of the stable 
region, since our estimate for i?*2 is only a lower limit. 
Thus, once transfer starts, it generally continues until 
both stars overflow their Roche lobes. 

3.5. The Fate of Unstable Mass Transfer Binaries 

What is the fate of a mass-transfer binary in which the 
primary cannot shrink inside its Roche lobe before the 
secondary also overflows its Roche lobe? Such a system 
will become a contact binary, or possibly even a com- 
mon envelope system if the semi-major axis decreases 
enough. Equation ([8]) implies that in a binary that trans- 
fers enough mass to reach q = 1, the semi- major axis 
shrinks by a factor of 2'*qQ/(l -I- go)'*- For go — 0.1, this 
is a factor of 10, and almost certainly produces common 
envelope evolution, or possibly even a merger. However, 
for qo = 0.5, the change in separation is only 20%, so a 
contact binary seems the more likely outcome. 

Once a contact binary or common envelope forms, over 
a thermal time the envelopes of stars will equilibrate to 
the same specific entropy. If the primary has a radia- 
tive core it will not equilibrate, but the radiative core 
may well disappear temporarily because as it expands 
and cools adiabatically, its opacity will increase, and it 
may transition back to convective instability. If this hap- 
pens, it will reach the same specific entropy as the two 
envelopes. The result will be an entirely isentropic con- 
tact binary. Since the two stars are isentropic, gas will be 
equally strongly bound to each of the two potential wells, 
and the system will approach equal masses. As the stars 
radiate they will continue contracting, their envelopes 
will shrink, and eventually contact will cease before the 
stars reach the ZAMS. The result will be a main sequence 
detached binary with a mass ratio close to unity. 

Note that this behavior is different than that com- 
monly observed in main sequence or post-main sequence 
contact binaries, which do not approach equal masses. In 
those systems, the two stars are chemically different (for 
example one may have consumed all of its hydrogen), or 
one or both may contain degenerate cores that will not 
reach the same specific entropy as the envelopes. It is 
these differences that produce unequal masses. 

After contact ceases and the stars contract to the 
ZAMS, the orbit will circularize and the binary will 
become synchronous. The characteris tic ti mescale for 
the fo rmer has been estimated by iGoodman fc DicksonI 
(|1998l ) (their eq. 15). Taking representative parameters, 
the circularization timescale for an equal mass binary 
with separation a ^ 0.1 AU, and constituent masses 
of 10 Mq and radii R - 0.02 AU, is icirc ~ 10^ yr 
(i?/0.02 AU)3/2[(a/i?)/5]2i/2. The very strong depen- 
dence on the ratio a/R makes this estimate quite un- 



certain: a factor of two change in a/R increases tdrc by 
a factor of 10^. Importantly, however, the timescale is 
shorter than or comparable to the characteristic thermal 
time. The synchronizatio n timescale (tsvn) can be est i- 
mated from the work of Goldrci ch fc NicholsonI (fl98l . 
Using the same parameters as for the estimation of ^circ, 
we find a synchronization timescale of tgyn ^circ with 
again a similar and very strong dependence on the mag- 
nitude of a/R. Although uncertain, these estimates in- 
dicate that close massive binaries shold be synchronous 
and on circular orbits relatively early in their main se- 
quence lifetimes. 

4. DISCUSSION 

A primary result of this work is that there is a critical 
a below which RLOF occurs (eq. [2l Figure [2]). Thus, for 
massive binaries that form at a separation of a few tenths 
of an AU, RLOF and mass transfer occur before either 
star reaches the ZAMS. For a relatively wide range of 
parameter space, mass transfer is unstable and we argue 
that the binary mass ratio will be driven toward equality 
on a dynamical timescale. For high system accretion 
rates, the ratio of the radiative core mass to the envelope 
mass of the primary protostar (^o) can be ^ 0.5, and in 
some cases the system can stabilize for some values of 
q < 1. Somewhat unequal mass binaries are the result 
when the stars reach the ZAMS (Figured]). 

4.1. "Twin" Binaries 

A particularly curious feature of the massive binary 
population is the high proportion of systems with mass 
ratios very clos e to unity. The mass ratio of WR20a i s 
q = 0.99 ± 0.05 (jBonanos et al.ll2004 iRauw et al.ll2005f) . 
and that of D33J013 346.2-H304439.9 is g = 0.90 ± 0.15 
(jBonanos e t al."2006'). Thirteen of the 21 detached bina- 
ries in the JHilditch ct al. (2005) sample have q > 0.85, 
and five have best fit values o f q consistent with g = 1.0 
(jPinsonneault &: Stanekl[2006f l. Due to its small siz e, the 
statistical significance of this sample is unclear (jLucvl 
120061 ). but a priori it seems quite unusual that almost 
25% of the eclipsing binaries in the SMC with orbital 
periods shorter than 5 days should have mass ratios con- 
sistent with 1.0. Nor does it seem likely a priori that 
the most massive binary known should have a mass ratio 
within 5% of unity. 

Ac cr etion from a circumbinary disk, as in the model 
iBatd (|2000D . provides a way of producing an anti- 
correlation between semi-major axis and binar y mass ra- 
tio, s uch as is observed even for low mass stars (jTokovininl 
l2000f ). However, to produce mass ratios near unity, this 
mechanism requires that binary protostars initially form 
with very small masses and then accrete many times 
their original mass from a circumbinary disk. While 
such low mass seeds may form in the cold, lower den- 
sity cores from which low mass stars for m, both analytic 
work ([Krumholz 2006) and simulations (jKrumholz et al.l 
|2007( ) suggest that radiation feedback largely prevents 
the formation of small "seed" st ars in massive protoste l- 
lar cores. In the simulations of iKrumholz et all (|2007| ). 
secondaries may form by fragmentation out of a mas- 
sive unstable disk (with subsequent migration), or via 
capture of a fragment formed elsewhere in the core, but 
in cither case the primary at the time of binary forma- 
tion is at least a few Mq in mass, and the secondary at 
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least a few tenths of Mq . Such a sy stem hkely could 
not reach a mass ratio of unity via the iBatd mechanism. 
However, as we have shown here, mass transfer between 
massive accreting protostars provides a natural, and for 
close binaries inevitable, mechanism for attaining mass 
ratios near uni ty. T hat said, the circumbinary accretion 
mechanism of iBatd probably operates in tandem with 
mass transfer — particularly at lower masses — bring- 
ing binaries closer together and pushing mass ratios to 
larger values, before mass transfer occurs and produces 
q of nearly 1.0 for tight binaries. 

4.2. Observational Predictions 

These calculations lead to definite predictions about 
the properties of massive binaries, which will be directly 
testable against future samples of binary systems larger 
than those available today. Since there is both a critical 
mass and a critical semi-major axis required for RLOF 
to occur, we predict that true twins, systems with mass 
ratios q > 0.95, should be significantly more common 
among stars with masses ^5 — 10 Mq with semi-major 

axes ^ 0.25 AU than among a binary population with ei- 
ther lower mass or larger separation. These twins should 
be in nearly zero eccentricity orbits, since mass transfer 
and evolution into a contact binary will circularize orbits 
rapidly. 

Furthermore, systems with small values of go that be- 
gin mass transfer will experience large reductions in the 
system semi-major axis. For example, a system with 
go = 0.1 will reduce its separation by a factor of 10, 
likely producing a merger. This may lead to a deficit of 
low mass companions to massive stars at separations sig- 
nificantly smaller than a tenth of an AU. Unfortunately, 
predictions regarding close binaries with small q are very 
hard to test observationally because such systems are 
difficult to detect. 

Massive stars should generally have larger mass ra- 
tios in clusters with high surface density, since these 
likely formed from higher pressure gas and thus pro- 
duced larger accretio n rates onto the stars within them 
(|McKee fc Tanll2003[ ). However, this effect is likely to be 
rather weak, since the critical semi-major axis only varies 
with the quarter power of accretion rate, and very high 
accretion rates make it easier for stars to remain unequal 
mass after the onset of transfer because they decrease the 
fraction of the star's mass that goes into the extended en- 
velope. As a result, there may be fewer true twins in very 
high surface density systems, even if there are more mas- 
sive stars with mass ratios ^ 0.5. In any event, due to 
weak dependence on the accretion rate, we consider this 
possibility less promising than searching for correlations 
of twin fraction with mass and semi-major axis. 

4.3. Primordial Stars 

One final note is that, although we have not discussed 
primordial stars, the mechanism we have discussed may 
well operate in them too. The critical ingredients for 
mass transfer are a phase of deuterium shell burning to 
produce large radii, and a close companion onto which 
to transfer mass. The binary properties of primordial 
stars are completely unknown, so we cannot comment 
on whether the second condition is likely to be met. 
However, a phase of deuterium shell burning does seem 



likely. Pr imordial stars proba bly form at high accre- 
tion rates (iTan fc McKed[200l . giving them large radii 
(jStahler et al.l 119861 )! In present-day stars, a radiative 
barrier forms because opacity in a stellar interior de- 
creases with temperature. The opacity source is primar- 
ily free-free transitions of electrons, and the availability 
of electrons for this process does not depend strongly on 
metallicity in an ionized stellar interior. Thus, primor- 
dial stars likely form radiative barriers much like present 
day ones, and undergo deuterium shell burning. 

5. SUMMARY 

Massive, rapidly accreting protostars can reach radii of 
tenths of an AU during their pre-main sequence evolu- 
tion, largely because they undergo a phase of deuterium 
shell burning that swells their radii. During this evolu- 
tionary phase, massive protostars with close companions 
will overflow their Roche lobes and transfer mass. Such 
transfer is always from the more massive, rapidly accret- 
ing star to the smaller one, since radius increases with 
both mass and accretion rate. Using simple protostellar 
structure and evolution models, we evaluate the range 
of separations and mass ratios for which mass transfer 
is expected to occur, and compute the likely outcome of 
the transfer. 

We find that, for the expected accretion rates in mass 
star- forming regions, binaries at separations of several 
tenths of an AU or closer will undergo mass transfer, with 
some dependence on the exact accretion rate and the 
initial mass ratio. This process always pushes binaries 
toward mass ratios of unity. For some accretion rates 
and initial mass ratios, mass transfer will either be stable 
initially, or will terminate on its own before reaching q = 
1. For many systems, though, it will only halt when the 
two stars form a contact binary. The stellar envelopes 
in such systems will rapidly reach almost equal masses 
and specific entropies, and then the stars contract onto 
the main sequence, forming a ZAMS binary with a mass 
ratio very close to unity. 

Pre-main sequence mass transfer represents a hereto- 
fore unknown phase of binary star formation and evo- 
lution, one that has likely affected a significant fraction 
of massive spectroscopic binaries. It provides a natu- 
ral explanation for the puzzling phenomenon of massive 
twins, high mass binaries with mass ratios that are con- 
sistent with q = 1.0. We predict based on this finding 
that twins should be significantly more common among 
stars ^ 5—10 Mq in mass at separations S 0.25 AU than 
among either less massive or more distant stars, and that 
mass ratios should generally increase with binary mass 
and decrease with separation. The weak dependence of 
Ocrit in equation ^ on the system mass accretion rate 
and the initial mass ratio suggests that this result should 
not depend significantly on formation environment. Sur- 
veys such as OGLE that detect statistically large samples 
of massive binaries are rapidly making these predictions 
testable. 
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